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Abstract. Recently helioseismic observations have revealed the pre- 
sence of a shear layer at the base of the convective zone related to the 
transition from differential rotation in the convection zone to almost uni- 
form rotation in the radiative interior, the tachocline. At present, this 
layer extends only over a few percent of the solar radius and no defini- 
tive explanations have been given for this thiness. Following Spiegel and 
Zahn (1992, Astron. Astrophys.), who invoke anisotropic turbulence to 
stop the spread of the tachocline deeper in the radiative zone as the Sun 
evolves, we give some justifications for their hypothesis by taking into 
account recent results on rotating shear instability (Richard and Zahn, 
1999, Astron. Astrophys.). We study the impact of the macroscopic 
motions present in this layer on the Sun's structure and evolution by in- 
troducing a macroscopic diffusivity Dt in updated solar models. We find 
that a time dependent treatment of the tachocline significantly improves 
the agreement between computed and observed surface chemical species, 
such as the ^Li and modify the internal structure of the Sun (Brun, Turck- 
Chieze and Zahn, 1999, Astrophys. J.). 

to appear in the Annals of the New York Academy of Sciences, Vol 
898. 



1. INTRODUCTION 

The presence of a shear layer connecting the differential rotation of the convec- 
tive zone to the solid rotation of the radiative zone is now well established by 
helioseismic inversions (Figure 1 and [1]). Both its location and width are more 
and more constrained and seem to be 0.691 it 0.004 Rq and less than 0.05 Rq, 
respectively. Today, there are several hydrodynamical or MHD descriptions of 
this shear layer and its extension but none is definitive [2], [3], [4]. In these 
models the motions in the tachocline are either turbulent or laminar, involve 
magnetic field or are purely hydrodynamical. In this paper we discuss the basic 
ideas supporting the dynamical description of this transition layer. First, in 
Section 2 we recall the physical processes acting in this layer and summarize the 
different approaches with an emphasis on Spiegel and Zahn's description invok- 
ing a nonlinear anisotropic turbulence. In Section 3, using the prescription of 
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Spiegel and Zahn [2] for the amplitude of the vertical velocity in the tachocline 
and Chaboyer and Zahn [5] for the chemical mixing and evolution, we build 
solar models including a macroscopic difFusivity Dt, which are compared to the 
most recent helioseismic data and surface abundance observations for ^Li, ^Be 
and ^He/^He ratio. Finally, we summarize our results and conclude in Section 
4. 



2. THE SOLAR TACHOCLINE 

The transition layer between the convective and radiative zones (Figure 1) plays 
a crucial role in our understanding of stars such as the Sun because it simulta- 
neously involves several physical features, such as: 

• the strong shear associated with the transition from differential to solid 
rotation, which may generate turbulence, 

• a turbulent interface with the convection zone above, which may produce 
internal waves [6] , [7] , 

• the possible presence of magnetic field, of fossile origin [3], or linked with 
the 11-year cycle [8], 

• the proximity of the thermonuclear burning zone of lithium 7 and beryllium 
9. 
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Figure 1. Solar rotation rate inferred from MDI data (aboard SOHO) 
as a function of radius at three latitudes, 0°, 30° and 60° (the vertical line 
represents the base of the convection zone) [9] . 

In this section we shall concentrate on some specific points concerning the 
tachocline and we will not deal with the full complexity of this transition layer. 
The presence of a latitudinal differential rotation at the base of the convection 
zone induces a latitudinal temperature gradient n{9) VT{9). Without any 
limitating process this temperature gradient will diffuse inwards, on a thermal 
diffusion time scale, enforcing differential rotation deep into the radiative interior 
of the Sun. As we already stated, however, helioseismic observations indicate a 
uniform rotation profile in the radiation zone, and thus we have to find which 
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processes could hinder this diffusion. 

There are different possibihties and first and foremost is the stable stratification 
of the radiation zone. This effect would indeed slow down the spread of the 
tachocline, but in spite of that, the layer would still extend to one third of the 
radius in the present Sun, as estimated by Spiegel and Zahn [2], which is far too 
much. Thus another process must be invoked to explain the observed thinness 
of the tachocline. 

One possibility is the presence of a magnetic field in the radiation zone, as advo- 
cated by Gough and Mclntyrc [3]. Their model is promising, but it has not yet 
been worked out in detail: it remains to be seen how the poloidal field threads 
into the convection zone, and avoids imposing differential rotation throughout 
the radiative interior. 

Another possibility has been suggested by Spiegel and Zahn [2]. If the latitudinal 
shear is unstable, it could generate anisotropic turbulence which would tend to 
reduce the differential rotation, and hence prevent the spread of the tachocline. 
It is not clear whether the solar tachocline is linearly unstable: when one applies 
the criterion derived by Watson [10] to an angular velocity profile of the type 
ri oc 1 — «! sin^ A, where A is the latitude, linear instability requires cti > 0.29, 
which is larger than the solar value ai ^ 0.25. It appears, however, that a law of 
the type fl cc 1 — ai sin^ A — 02 sin^ A, closer to the latitudinal dependence drawn 
from helioseismic inversions (see Fig. 1) would be more sensitive to such insta- 
bility [11]. Furthermore, a toroidal magnetic field of sufficient strength would 
also act to destabilize the flow, as shown by Oilman and Fox [4]. In any case, the 
Reynolds number is so high that such a differential rotation would be liable to 
finite amplitude instability, as can be infered from laboratory experiments [12]. 
It has been argued by Oough and Mclntyre [3] that turbulence does not nec- 
essarily reduce the shear of differential rotation: in the Earth's atmosphere the 
transport of angular momentum would even imply a negative viscosity. But 
there the transport is achieved mainly through Rossby waves, and it is not clear 
whether this can be applied to the solar tachocline. 

Admittedly, this important issue is still a matter of debate, and it may be set- 
tled only by comparing the models' predictions with the observed properties of 
the Sun. This is why we have chosen to draw all observable consequences from 
the model which has been worked out in sufficient detail to allow such a test, 
namely the turbulent tachocline proposed by Spiegel and Zahn [2]. We will show 
that the mixing induced in this model improves both the sound speed profile 
(reduction of the peak below the convection zone) and the surface light element 
abundances, confirming the need of introducing macroscopic processes in solar 
models [13]- [15]. 

3. MIXING IN THE SOLAR TACHOCLINE: PHYSICAL DESCRIP- 
TION 

Macroscopic mixing may be treated in solar models by adding an effective diffu- 
sivity Dt in the equation for the time evolution of the concentration of chemical 
species. To establish this coefficient for the tachocline, we use the description by 
Spiegel and Zahn [2], where anisotropic turbulence is responsible for stopping 
the spread of the layer. This anisotropic diffusion will also interfere with the 
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advective transport of chemicals. Chaboyer and Zahn [5] have shown that the 
result is a diffusive transport in the vertical direction. Using their result, Brun, 
Turck-Chieze and Zahn [16] derived the following expression for the macroscopic 
diffusivity: 

Dxir) = j^VH {-—\ Ol exp(-2C) cos^(C) + higher order terms (1) 

4Uo \fhcz J 

with Q/i = 04/0, characterizing the differential rotation rate, ^4 = 4.933, 
C = lJ-4{rhcz — r)/d a non-dimensional depth, 

d = n,A2^/Nf/\AK/uHf'^ (2) 

a length related to the tachocline thickness h (e.g h ~ c?/2), r^cz the radius, 17 
the angular velocity and K = xl 9^%, the radiative diffusivity at the base of the 
convective zone. The horizontal component of the macroscopic diffusivity is 
assumed to be equal to the horizontal viscosity vu. In our solar models, we treat 
h (hence d) as an adjustable parameter, chosen to agree with the helioseismic 
determination of the tachocline thickness h < 0.05Rq [1]. With the latitudinal 
dependence of the angular velocity at the base of the convection zone deduced 
from Thompson et al. [17], Qhcz/'^'^ = 456 — 72x^ — 42x^ nHz (with x = sin A), we 
have reestimated the coefhcient = — 1.707 x 10~^, as well as the ratio between 
the rotation in the deep radiative zone and the equatorial rate Q/^Iq = 0.9104. 
The prediction by Gough and Sekii [18] for the latter, who consider instead 
the magnetic stresses, is ~ 0.96; presently, the seismic observations suggest 
a rotational ratio of 0.94 it 0.01 [1], which is intermediate between these two 
theoretical estimates. 

An analysis of the dependence of our d and Dt with the global and differ- 
ential rotation rates yields 

Dt(xuh( — ) ocnulj\n/nf, docn^/^/u]{\ (3) 

\ ^bzc / 

where we have used equations (|l|) and (^) . Assuming that the turbulent viscosity 
is proportional to the differential rotation (i.e., uh oc (l), as suggested by the 
laboratory experiments [12], and introducing the dependence of the differential 
rotation on rotation observed by Donahue, Saar and Baliunas [19] {Cl oc Q,^'"^^^'^), 
we finally obtain the following scalings 

oc d OC (4) 

We conclude that the tachocline mixing was stronger in the past both because 
that layer was thicker and because the diffusivity was larger. We render the 
mixing in the tachocline time dependent, through Dxi^it)) and d{0,{t)), by 
using the spin-down law 17 oc t^^/^ which was deduced by Skumanich [20] from 
the rotation rate of stellar clusters of different ages. 
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4. RESULTS 



Starting from the reference model of Brun, Turck-Chieze and Morel [21] built 
with the CESAM code [22], we introduce for this study the nuclear reaction 
^Li(p,a)^He proposed by Engstler et al. [23] and the coefficient Dt (Eq. 1) in 
the diffusion equation of chemical species, and we follow the time evolution of 
the solar model from the pre-main sequence (PMS) until 4.6 Gyr. The results 
are shown in the Table 1 and Figures 2-5 (see Ref. [15] for a more detailed 
discussion). We use a tachocline thickness h of 0.05 or 0.025 Rq, N of 100 or 
25 /iHz and Q, of 0.415 /iHz. Our standard model [21] has a surface abundance 
for helium of 0.2427 in mass, corresponding to an He difffusion of 10.8%. This 
value of Yg is a bit too low if we compare with the Basu and Antia [24] value 
for the OPAL equation of state [25], = 0.249 ± 0.003. 

Table 1. Surface abundance variation of ^He/'^He during the last 3 Gyr, 
surface abundances of *He and heavy elements Z, and abundance ratio ini- 
tial/surface for ^Li and ^Be from observations and for solar models at the 
solar age" 
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Ct 




< 0.1 




0.1 


0.1 


0.1 


0.1 


0.1 


0.05 


N (/iHz) 






100 


25 


100 


25 


25 


25 


(3He/4He), 


max 10% 


2.28% 


2.14% 


2.01% 


2.11% 


2.0% 


2.02% 


2.07% 




0.249±0.003 


0.2427 


0.2452 


0.2473 


0.2455 


0.2477 


0.2509 


0.2464 


(Z/X). 


0.0245±0.002 


0.0245 


0.0245 


0.0245 


0.0245 


0.0245 


0.0255 


0.0245 


^Lio/^Li, 


~100 


~6 


~8 


~22 


~12 


~91 


~134 


~89 


9Bco/9Bcs 


1.10±0.03 


1.115 


1.093 


1.086 


1.093 


1.118 


1.125 


1.092 



"Subscripts t for time dependent models, tz for time dependent model with Zq =2^'^ = 0.01959 
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Figure 2. Radial profile of the difference of ^He composition between the 
initial and present values for the reference solar model including only mi- 
croscopic diffusion {solid line) and solar models where we add a macroscopic 
mixing due to the presence of the tachocline: coefficient A (dash) and B {dash 
dot) (see Table 1 for the model characteristics). 

When introducing our diffusive coefficient, we mix helium back into the 
convection zone, inhibiting the microscopic diffusion up to 25% and producing 
a photospheric ^He^ = 0.2473 (cf. Table 1, models A and B). As expected, 
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Figure 3. Squared sound speed difference between GOLF+MDI data and 
ttie reference model (solid line) or models including a turbulent term: coef At 
{d = 0.1, N = 100) (dash), Bt {d = 0.1, N = 25) (dash dot), Ct {d = 0.05, 
N = 25) (dash three dots), and one model with Bt and Zq = Z'^q^ = 0.01959 
(i.e., Btz) [long dash). 



the composition profile is smootlier and flattens over the distance h below the 
convective zone (see Figure 2). 

The effect on the sound speed is displayed in Figure 3. When the macro- 
scopic transport is neglected, the squared sound speed difference reveals a peak 
just below the convection zone, coinciding with the tachocline (solid line). Macro- 
scopic diffusion acts to reduce this peak, but when one recalibrates the model to 
yield, e.g., the present abundance of heavy elements (Z/X = 0.0245=b0.002), the 
effect is rather minor. On the other hand, if the heavy elements are let free to 
adjust, within the observational uncertainties, the peak is completley removed, 
leaving only a broad bump culminating at 0.6 Rq, which presumably is due to 
another cause (long dashed line). 
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Figure 4. ^Be and ''Li radial profile (approximately superimposed the tem- 
perature scale) for several models: respectively, reference (solid line), coeffi- 
cient At (dash), coefficient Bt (dash dot) and coefficient Ct (dash three dots). 

The two light elements ^Li and ^Be are extremely sensitive to mixing pro- 
cesses occuring in stars because their nuclear burning temperatures are rather 
low (respectively, 2.5 10*^, and 3.2 10^ K) [26]. The new observational constraints 
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can only be satisfied if those chemical species are mixed in a rather thin layer 
below the convective zone, in order to preserve ^Be, which is very little depleted 
according to Balachandran and Bell [27]. This is the case with our tachocline 
model. However, if the mixing had proceeded in the past at the same rate as 
in the present Sun, ^Li would have been depleted only by a factor ~ 4, which 
is insufficient to account for the photospheric lithium abundance ([28] and ref- 
erences therein). But the thickness of the tachocline and the strength of mixing 
have been larger in the past, when the Sim was rotating faster. This effect is 
included in the models labeled with the index t (as Bt), whose evolution was 
calculated with the time-dependent diffusivity of Eq. (4) . 

3.51 J 
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Figure 5. Time dependent depletion of ^Li for several solar models: no 
microscopic diffusion (dash), with microscopic diffusion (solid line) and with 
mixing in the tachocline thickness: coefficient A {dash dot), B {dash three 
dots), then with time dependent mixing At {dots), Bt {long dash), Ct {thick 
dash dot) and Btz (Zq =T^q^ = 0.01959) {thick dash). We superimposed on 
the theoretical curves the open cluster observations (adapted from Vauclair 
and Richard [14] and cluster age uncertainties deduced from Lebreton et al. 
[29]). 

In Table 1 we give the initial over present ratio of ^Li and ^Be and show in 
Figure 4 the radial profile of ^Li and ^Be normalized to the surface abundance. 
We clearly see that the mixing process modifies the distribution of lithium but 
not that of beryllium (exception being the flat plateau for the mixed models in 
comparison with the "pure" diffusive one). With the coefficients B more ^Li IS 
burned than with and we also sec that the time dependence (models with 
index t) improves the ^He surface abundance as well as the ^Li depletion, where 
a value of ~ 100 is obtained without destroying ^Be or increasing too much the 
^He/^He surface ratio over the past 3 Gyr, as deduced by Geiss and Gloeckler 
[30] from meteorites and solar wind abundance measurements (see Table 1). 

In Figure 5 we show the lithium depletion occurring during the Sun's evo- 
lution for the different models presented, plus a model without any diffusion. 
Clearly, only the diffusive models including mixing in the tachocline yield a sub- 
stantial depletion during main sequence evolution, in agreement with the obser- 
vations (superimposed with their inherent dispersion on the theoretical curves). 
Note that the strong time dependent mixing with = 25 (models Bt, Ct and 
Btz) presents a reasonable value of the solar ''Li depletion (~ 100). 
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However the lithium depletion during the PMS is probably overestimated 
due to the crude spin-down law we have adopted. A more detailed analysis of 
these phases is under study, including metallicity effects and more appropriate 
angular momentum evolution during this phase. 

Our results show the interest to follow together the photospheric abundance 
of the four elements ^He, ''He, ''Li, '^Bc, and to examine their sensitivity to the 
microscopic, as well as the macroscopic, processes. This study encourages the 
introduction of macroscopic processes in stellar evolution models, and demon- 
strates the crucial role of the thin tachocline layer below the convective zone. 
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